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Electronic Characteristics of 
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Na x HfNCl (T c =25 K) 
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Local density functional results are presented for the 
electron-doped metallochloronitrides yl^ZrNCl and AajHfNCL 
A=Li or Na, which superconduct up to 25K. The alkali non- 
stoichiometry is treated in a virtual crystal approximation. 
The electronic structure is strongly two dimensional, espe- 
cially in the conduction band region occupied by the carri- 
ers, because the states are formed from the in-plane orbitals 
d xy ,d x 2_ y 2 of the metal ion and the Px,p y orbitals of the N 
ion. We predict a change of behavior at a doping level of 
z=0.3. 



INTRODUCTION 

High temperature superconductivity (HTS) in the lay- 
ered cuprates continues to puzzle even after a dozen years 
of intense scrutiny. Even outside of the class of HTSs, 
however, the appearance of high T c (in the pre-cuprates 
sense) also is baffling. The highest T c is achieved in the 
fullerides, with T c =40 K reported, Q which arc three 
dimensional (3D) materials with, however, a great deal 
of zero dimensional (cluster) character. T c ~ 35 K is 
achieved in Bai-^K^BiOa, also a 3D system, with only 
some flat portions of Fermi surface to bring in questions 
of low dimensionality. There is not yet any clear theo- 
retical accounting for the magnitude of T c in these two 
systems |^|,|| although some phonon modes in the ful- 
lerides are strongly coupled to the conduction electrons. 
The A15 superconductors (T c up to 23 K) present a class 
in which the relatively high value of T c is understood in 
terms of strong electron-phonon coupling and a high den- 
sity of states (DOS) [N(E F ] at the Fermi level E F . §j 

A new, 2D material has recently been added to the 
list of (non-HTS) superconductors with unusually high 
T c . Yamanaka et al. || discovered superconductivity up 
to 12 K in Li-doped ZrNCl, and up to 25.5 K in Na- 
doped HfNCl. Shamoto et al. j7|,|| have reproduced 
the superconductivity and determined the crystal struc- 
ture of superconducitng materials using synchrotron X- 
ray diffraction data. We provide here the first calculation 
of the electronic structure of this system, considering in 
particular variation with doping level. 



DESCRIPTION OF CALCULATIONS 
Electronic Structure Methods 

We have applied the local density approximation || us- 
ing two methods: the linearized augmented plane wave 
(LAPW) method for the bands and density of states 
(DOS) that we show, and the plane wave ultrasoft pseu- 
dopotential method jll] for the phonon frequencies. Both 
methods produced essentially the same band structures. 

To account for the doping we have used a virtual crys- 
tal treatment. For example, for x fractional occupancy 
of a site by Na (ten core electrons and one valence elec- 
tron), that site is fully occupied by a nucleus of charge 
10+ir with the corresponding number of electrons. As ex- 
pected, Na is ionized with the electrons going into bands 
that are primarily metal d bands, so the virtual crystal 
treatment should be reasonable. We will find that it is 
important to treat this charge self-consistently, because 
non-rigid-band behavior appears not far from the region 
of immediate interest . 



Structure 

We have used the structures of Shamoto et al. [fj],[| for 
these materials. Because a Cl-HfN-HfN-Cl slab is only 
weakly coupled to its neighboring slabs, we expect that 
the stacking sequence (for example, whether hexagonal 
or rhombohedral) is not important, and indeed that was 
found to be the case. One such structure is shown in Fig. 
1, and the slab structure is the same for Hf and Zr, with 
only the stacking sequence varying with doping. |12| 



DISCUSSION OF RESULTS 

The band structures are strongly two dimensional, so 
we discuss only the in-plane dispersion here. The un- 
doped compounds have a band gap of 1.7-1.8 eV, be- 
tween valence bands that are 2p states of CI and N, and 
conduction states that are mostly Hf or Zr. Thus the 
formal ionic description (Hf,Zn) 4+ N 3 ~Cl~ is reasonable, 
but obscures the strong metal-N hybridization. 





FIG. 3. Atom (left panel) and orbital (right panel) pro- 
jected DOS for Nao.25HfNCl. In the right panel only the 
conduction bands are pictured. 



FIG. 1. Rhombohedral structure of Li^ZrNCl determined 
by Shamoto et al. [7,8] The stacking along c can change with 
doping, but the basic slab structure structure is unchanged 
with doping. Li sites are only fractionally occupied 
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FIG. 2. Plot along symmetry lines of the band structure 
of Nao. 25 HfNCl. Points are F = (0,0,0), K = (2/3,1/3,0), 
M = (1/2,0,0), A = (0,0,1/2), in units of the hexagonal 
reciprocal lattice vectors. The lack of dispersion along T-A of 
the lower conduction bands reflects the strong 2D character 
of the important conduction band. 



In Fig. 2 the band structure of Nao.2sHfNCl is shown 
along the hexagonal symmetry directions. Only a single 
band is occupied by the doped carriers, and the band 
structure is very similar to that of the undoped com- 
pound, i.e. a rigid band picture is good. The projected 
DOS of Fig. 3 illustrates several things: (1) the occupied 
band is a hybridized planar band of Hf d xy , d x 2_ y 2 and N 
p x ,p y character; (2) Hf d xz , d yz and N p x character do not 
appear until the band flattening 0.3 eV above Ep in Fig. 
2; (3) the Hf d z 2 character is higher still. This behavior 
is completely different from what Woodward and Vogt 
p3[ calculated for a bilayer structure reported by Juza 
and collaborators, Jl4| in which two Hf ions are bonded 
across the bilayer, rather than the Hf-N bonding across 
the bilayer as in the structure determined by Shamoto. 
Thus electrons are not doped into a very flat band, which 
almost certainly would be strongly correlated, but rather 
into a broad band with a light mass m* rs 0.6m. 

The other material parameters we obtain (for x=0.25) 
are 

• Fermi velocity yp = 3 x 10 7 cm/s 

• Drude plasma energy ft p — [Aire 2 N(Ep)v F J 1,/2 = 1.5 
eV 

• Gap 2A « 3.5 k B T c « 7.5 meV 

• Coherence length £ = ^ » 180 A 

• London penetration depth A = jf-= 1300 A 

• k = A/£ = 7^ Type II superconductor. 



FIG. 4. Plot of the distorted circular Fermi surface of 
Nao.25HfNCl shown in two hexagonal Brillouin zones. The 
surface is centered at the zone corner point K; since there are 
two such inequivalent points, there are two Fermi surfaces. 

In Fig. 4 we present the Fermi surface for 2=0.25 
doping level. It consists of threefold distorted circles 
at the zone corner (K) points. There are two inequiv- 
alent such points, each with mean radius kp — |f • On- 
Fermi-surface scattering processes will be dominated by 
(1) small Q< 2fcp intraband scattering, and (2) large 
Q^K interband scattering. [Note that the vector con- 
necting two neighboring zone corners K is also the vector 
K.] 

We have calculated the frequencies of the three fully 
symmetric (Ai s ) Raman active vibrational modes of 
pristine ZrNCl and HfNCl at the zone center, cor- 
responding to modulation of the internal coordinates 
z Zr/Hfi z n, zci- The frequencies are 586, 334, and 202 
cm -1 and 604, 368, and 191 cm -1 respectively for the 
two compounds. The highest frequency is almost pure N 
motion. 



DOPING 

As these systems are studied further, the question of 
the effect of doping will be very important. We have 
performed virtual crystal calculations for several doping 
levels, with results for x = 0.25, 0.35 and 0.45 shown in 
Figure 5. There are continuous rigid-band-like changes 
up to just above x = 0.30, above which bands with dif- 
ferent character (specifically, out of plane orbitals) begin 
to become occupied. At this point in doping the proper- 
ties should show a change in slope (plotted versus x). At 
cc=0.35 Ep coincides with a band that is very flat band 
along T-K and has a small mass along T-M as well, i.e. 
a peak in the density of states. This range of doping 
should prove very interesting, and whether T c rises or 
falls will reveal important characteristics of the pairing 
mechanism. 



FIG. 5. Band structure plots of Na* HfNCl, for z=0.25, 
0.35, and 0.45 (left to right). Around £=0.35 the Fermi 
level enters a region where the bands have different charac- 
ter, and already the conduction band dispersion in visibly 
non-rigid-like. At £=0.45 there is a large Fermi surface cen- 
tered at F as well as the ones at K. 



SUMMARY AND ACKNOWLEDGMENTS 

It is surprising that such low mass, low DOS materials 
such as those described here can superconduct up to 25 
K. Since there is no indication, experimentally or theoret- 
ically, of strong correlation effects in the range of doping 
reported by the superconductors, a possible candidate 
for pairing mechanism could be phononic. Given the 2D 
character and the strong nesting that might introduce 
characteristics of ID behavior, purely electronic pairing 
mechanisms should also be considered. Given the appar- 
ent extreme two dimensionality of the crystal, it will be 
important to establish that coherent supercurrent actu- 
ally will flow perpendicular to the layers. 
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No. N00014-97-1-0956 and National Science foundation 
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